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ABSTRACT 

This module is part of a series designed to be used 
by life science students for instruction in the application of 
physical theory to ecosystem operation. Most modules contain computer 
programs which are built around, a particular application of a 
physical process. This module is used to introduce the biology 
student to differential calculus^ ^ br?inch of mathematics which is 
being used more often in recent ecological and physiological models. 
Since biological systems are dynamic, their mathematical models must 
describe rates of change of relevant variables , a process which 
requires calculus. The module introduces modeling in biology by 
reviewing differential calculus and using only examples from life 
sciences- A problem set reviews information in the text proper and 
presents additional information not included in the text- An 
associated computer program, DIFF, uses graphics to check the user's 
own calculations and to demonstrate the validity of general 
solutions. (Author/CS) 
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PREFACE 

* Calculus has quietly invaded several areas of biology in the 

last few years, reflecting a greater desire for precise explanations 
of biological phenomena. Mathematical modeling of rates of dynamic 
processes requires calculus. This module introduces modeling in 
biology by reviewing differential calculus using only examples from 
life sciences. The problem set should be worked since several ideas 
are presented which are not in the text proper. An associated com- 
puter program, DIFF, uses graphics to check the user's own calculations 
and demonstrate the validity of general solutions. Previous exposure 
to calculus is recommended but not required. 
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IlilTRODUCTION 

Differential calculus is being used more often in recent ecological 
and physiological models as data become more precise and the processes 
become better understood. Since most biological systems are dynamic, 
their mathematical models must describe rates of change, not just current 
values, of the relevant variables. Although most models consider changes 
over time, the techniques of calculus depend only on the mathematical 
function involved and thus any independent variable may be used, such as 
spatial dimensions, organism weight, temperature, etc. As a result, the 
implications from a model of one process may be applied to the model of a 
different process as long as the mathematical functions involved are the 
same. 

FUNCTIONS OF ONE VARIABLE 
Rates of Change 

The simplest graph of a dynamic relationship is a straight line. The 
equation for a straight line is 

y = mx + b (1) 

where y and x are variables, m and b are constants. An example of this 
relation is the oxygen uptake by the lobster. The oxygen consumption (y) 
depends on the oxygen concentration (x) in the surrounding environment, so 
that y is a function of x. A typical graph of this function in shown in 
Figure 1. 
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The equation for this function is 

y(x) = .004x + .016 
The number .004 represents the slope of the line; that is, the ratio of 
the change in y to the change in x. When x changes by 1 unit, y changes 
by .004 units. In this particular application, when the water O2 concen- 
tration increases by one ml/l, the lobster O2 consumption increases by 
.004 ml per hour^-gm body weight. The slope (m in Equation (1)) then 
represents the rate of change of y with x. 

When the graph is not a straight line, the function it represents is 

I 

more complex than above and the rate of change cannot be expressed so 
easily. Note that for each unit change in x in Fig. 1, y changes by .004, 
regardless of t^ie value of x. The rate of change is then constant . In the 
graph of Fig. 2, the rate of change is not constant. To see this, approxi- 
mate Fig. 2 by two connected tangent lines (Fig. 3a) and note that the 
slope differs with each line. As the approximation improves (Fig. 3b), it 
uses more lines and thus presents more slopes. Using an infinite number 
of lines, we would duplicate the curve (in Fig. 2) and have a slope that 
changes with each value of x. The slope then depends on x and clearly is 
not constant. In fact, one definition of the slope of a curve at a point 
is the slope of the tangent line at that point. 

Note that the slope is ambiguous at the points where two straight 
lines meet (Figure 3a), We say the slope is "undefined" at such "corner" 
points . 

In general, the rate of change of a function is also a function of x 

dy 

and possesses its own equation. The rate of change is denoted to 
reflect the ratio of the change in y to the change in x. When the graph 



X 



Figure 2. Example of a function with a changing slope. 





a. b. 

Figure 3, Approximations to curve in Figure 2 . 

is a straight line, the equation for y is 

y(x) = mx + b 
and the rate of change is 



dx 



= m 



The rate of change is called the derivative . Its functional form 
depends on the equation for y. For example, an empirical relation 
between oxygen consumption (Q) and body weight (W) is 

Q = SW^ 
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The derivative of this function is 
^ = (3)(2W) - 6W 

The graph of Q - 3W^ is given in Fig. 4a. The derivative at W « 1 represents 
the slope of the line tangent to the curve at the point* where W = 1, as 
shown in Figure 4b. At W = 1, the slope is calculated to be 

§.6(1,-6 





slope=6 



W 



0 1 2 3 W 



a. b* 
Figure 4. Deri vative as the slope of the tangent lin e* 

The number "3" in the formula Q = 3W^, as well as the "2" in the 
exponent, are empirically determined and differ with body size and species. 
The general formula is 



Q = aW 



The derivative of this power function is 



(2) 



-77^ = a»b»W 
dw 
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Table 1 gives the more common functions and their derivatives, 



Table 1. Derivatives of 


Elementary Functions 


y(x) 


dy/dx 


n 

X 


n-1 

nx 


X 

e 


X 

e 


Jin X 


1/x 


sin X 


cos X 


cos X 


-sin X 



For a more complete table, see any calculus text, or » 
any math handbook (see Bibliography) . 

Composite Functions 
When a function is composed of several simple functions its derivative 
can be evaluated in stages. In the simple cases where y equals the sum or 
product of two functions, f(x), g(x), the rules for differentiation 
(finding the derivative) are: 

y(x) =f(x) -Hg(x), S = f + ^ 

y(x) =f(x) g(x), g= f(x)f +g(x)f . 
For example, if y(x) = x^(x-l) + 2x, then 

d2=,2l(x-l) |(x^)^2 
dx dx dx 

= x^d) + (x-l)(2x) + 2 

When y is a function of a function, the chain rule provides the differen- 
tiation method: 

10 
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dy dy du 
dx du dx 



2x 

For the composite exponential function y « e , we have 



u(x) = 2x, y(u) = e" 



du * dx 



The refore 



dx 



The Gompertz growth curve is used occasionally to describe the popula- 
tion size (N) of some species as a fimction of time (t), and is given by the 
equation 

, -kt 

N(t) = ae" 

where a,b,k are constants. The derivative dN/dt then represents the rate of 
growth of the population. Here the chain rule is applied twice: 

N(t) = ae''^^^ u(t) = -be"^^ 



dN u 
^ = ae 
du 



Write u(t) as u = -be^^*"^ where v = -kt . Then 



dv , 'du du dv . , . 

dt du dv dt ^ 



Higher Derivatives 
The derivative dy/dx of a function y(x) is called the first deriva - 
tive of y(x). If we write 



then differentiating g(x) produces 

s - 

I 

which is called the second derivative of y(x), and is written 
dx^ 

Other notation used for the first derivative includes y* and y, for the 
second derivative, y** and y. Since the second derivative is also a function, 
it too can be differentiated to give the third derivative, and so on. The 
n^^ derivative is written (there is no general dot notation) 

^ n * ^ 
dx 

Critical Points 

The first and second derivatives can be used to determine three special 

points on the graph of the function, namely, the relative maxima, relative 

'minima and the inflection points. The relative maximum is easily visualized: 

the curve rises, reaches a peak, and then falls. The peak is the relative 

maximum. It is relative because the curve may rise even higher in a 
y^ 

different place on the graph. Similarly, the relative minimum constitutes 
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a low point on the curve. An inflection point is best illustrated by an 
example. 

The logistic growth function describing population size is 

N = N , - N (l+b)(l+be"'^^)""^ 

1+be 

where k is a growth coefficient, is the population size at t = 0 and N^(l+b) 
represents the carrying capacity of the environment. The growth rate is then 

^ = N (1+b) (-1) (1+be""''^)--^ (be"^^) (-k) 

G t O 

= N^(l+b)bke"^^(l+be"^^)"^ 

The coefficient k is always positive. In this example, we restrict b to be 
greater than 1. 

At a relative maximum, the tangent line is horizontal so the slope is zero. 
Then the maximum growth rate occurs when the derivative of the growth rate 
equals zero. 

N 



,slope at (t^,N^) 



dt 



or equivalently. 



= 0 



dfN 
dt2 



= 0 




inflection point (t^,N^) 



Figure 5. Maximum growth rate at the inflection point 

The second derivative of a functior equals zero at the point of inflection, 
where the curvature changes from cowing upward to ( irving downward , 

or vice versa. T?jen the maximum growth rate occurs when the population function 
N(t) is at its inflection point (see Fig. 5). Since the slope is decreasing 
(leveling off) following the inflection point, and increasing before the inflec- 
tion point, it is certainly maximal (steepest) ait that point. This point can 

13 
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also be viewed as the relative maximuTn on the graph of the growth rate, 
dN/dt (see Fig. 6). 

dN/dt 

zero slope 




Figure 6. Growth rate as a function of time . 
We now locate this point of maximum growth: 



0 = 



d^N 



at (t^, N^) 



0 = ■|^[N^(l+b)bke~^*'(l+be"^*')~^] at (t^,N^) 



-kt, -kt. 



-kt. 



0 = (l+b)bk e ■'•(be ■'•-Dd+be . 

-kt^ 

Since all factors are positive except (be -1) , then 

-kt 

0 = be -1 



Solving for t^ gives 



and then substituting into the original expression for N, 

-k(i£nb) 

N = N (l+b)(l+be ^ ) = N (l+b)/2 . 
X o o 

this says that the growth rate is highest when the population is one-half 
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of the carrying capacity. 

FUNCTIONS OF SEVERAL VARIABLES 
Partial Derivatives 

The models treated thus far involve functions of one variable. Oxygen 
uptake is given as a function of just the surrounding oxygen concentration. 
The population size depends only on time. A more complicated model, however, 
may involve many variables. Growth certainly depends on available food sup- 
ply in addition to time. Oxygen consumption also depends on more factors 
than ambient oxygen concentration. One such model is discussed by Bayne^ 
Thompson and Widdows (1973). 

The model begins with the equation 

§ = aW^ (3) 
at 

where C is the amount of oxygen consumed up to time t, and a, b and W are 
constants.* The notation dC/dt demands that we may be able to consider C 
only as a function of t. This is not always the case. Bayne, et al., 
studied mussels (Mytilus) with regard to the effects of food and temperature 
on oxygen consumption. One of their data sets gives values for the co- 
efficients a and b for winter vs. summer at two activity levels: 

Table 2. Oxygen consumption for Mytilus edulis 

Activity 

Parameter Season standard routine 

^ Winter 1.76 2.6A 

Summer 1.87 2.b4 

b Winter 0.724 0.774 

Summer 0.670 0.702 



^ dC 
Comparing (3) with (2) gives "J^ Q- 
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Standard activity represents the resting state. Routine refers to the post- 
feeding time period where some filtration (i.e. muscle activity) is occurring. 
From table 2, we recognize significant dependence of "a" on the level of 
activity and dependence of "b" on both activity and season. Let season and 
activity be denoted s and m, respectively. Then, in place of "a" and "b'' we 
write a(m) and b(s,m) to show the dependence on the variables m and s. Now 
C, a and b are dependent variables and W, s, m and t are the independent 
'Variables. We must now write (3) as 



The derivative notation is different from that in (3) to indicate more than 
one independent variable. This derivative is called a partial derivative 
and represents the rate of change of C with time while all other variables 
are held constant. Note that by holding all of the independent variables 
(except t) constant, we also hold a and b constant. So this partial deriva- 
tive is obtained by dif fe ren t i at in g the fun c t ion C wi th respe c t to t an d 
treating all of the remaining independent variables as constants* 
As a simple example, consider the function 



rb(s,Tn) 



(4) 




Then 




(t held constant) 



and 



x*2t (x held constant) 



1 
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A partial derivative is by nature merely one simple relation (of 

many) extracted from a complicated function. In (5), when both x and t vary, 

the graph of y vs . x vs. t is a ^hree-dimensional surface (figure 7). When 

2 

t is held constant, the graph (y vs. x) is a straight line (with slope t ); 
when X is held constant, the graph (y vs. t) is a parabola, as shown in 
Figures 8 and 9. These latter two graphs are much simpler than the surface 
of figure 7. Note that the straight line (y vs. x) is the far edge of figure 
7a, and the parabola (y vs. t) is the near edge of figure 7c. 




0 X 0 t 



Figure 8. dy/^x = t^ . Figure 9. ay/3t = 2xt . 

A more complicated example is the complete expression for (4) . The level 
of activity is based on the fraction of the maximal filtration rate. Then 
m=0 represents the ''standard" state, m=l gives the "active" state, and 
m=.4 is the "routine" state . Consider oxygen consumption during summer. 
Then 

I 

a(m) = 1.87(tiH-1) 
and (3) becomes 

1^ = 1.87(m+l)W^-'^^ (e>) 

o t 

Changes in "b" are not significant so that an average value, .7, can be 
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used. Since (6) involves four variables, and thus cannot be plotted, we 
revert to the notation of (2), i.e*, Q « • Then 

Q ^ 1.87(in+l)W^*^^ (7) 

This last expression is similar in form to (5) and its graph has a similar 
shape. Problem 3 discusses (7) in more detail. 

Critical Points in Three Dimensions 

Tlie extension of a critical point to functions of two variables is quite 

natural. A point P = (x ,y z ) is a critical point for the function 

o -^o, o 

2 = f(x,y) 

if 

M.ii. 0 

9x 9y 

at the point P. The classification of the critical point is, however, more 
complicated. Since there are three second derivatives, many cases could be 
considered: 

ill - il9f/3x) d^f ^ Hdf/dy) 

3x2 3x 3y2 3y 



3^f = d(^df/dy) _ ^3f/3x) = 3^f 
3xoy 3x " 3y " 3y3x 



This last "mixed*' second derivative can be evaluated in either order only 
if the function f (x,y) is continuous in x and y. The functions used in 
the examples wh ;ch follow are continuous so that the order of differentia- 
tion is arbitrary. Rather than considering all combinations of sign (+,0,-) 
in the second derivatives, we treat only three, which classify a relative 
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maximum, relative minimum, and a saddle point. Define 




1 

where 

P is a relative maximum if L > 0 and f (x ,y ) < 0. 

XX o o 

P is a relative minimum if L > 0 and f (x ,y ) > 0. 

XX o o 

P is a saddle point if L < 0. 
When L=0, the situation is ''undetermined" since its resolution is beyond 
the scope of this module. 

An interesting example is the function 

2 = x^ + y^ - 3xy + 15 

which seems to describe some of the properties of water falling across a 
rock face (Clow and Urquhart, 1974). Some of these properties are well 
known. The water will often dig potholes in the rock, especially if it 
falls onto a ledge. The corners and edges of the ledge eventually become 
rounded. We would then use a function which drops rapidly, levels off 
then drops steeply again. Figure 10 shows a three-dimensional computer 
plot of this function, looking across the origin into the positive octant 
(x>0, y>0, 2>0). The required shape is evident, with the pothole just 
beginning to form. In fact, the function does possess a relative minimum 
at the point (x, y, z) = (1, 1, 14). This model is discussed further in 
problem 6. 

The saddle point in the waterfall model is located at the point 
(0, 0, 15). The region around the saddle point represents the front 
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part of the ledge. It is displayed in the computer-drawn graph of figure 
11, expanded vertically to highlight the saddle shape. Some basic features 
of derivatives are shown here: 

a) The slope changes from point to point, 

b) The slope depends on the orientation of the tangent line. 
Thus, at a given point, the slope found by Sz/3x may be 
different from the slope found using dz/dy^ 

c) The slope at a relative maximum or relative minimum is 
zero , i.e. horizontal. 




Figure 10, Waterfall function . Figure 11. Saddle point of waterfall 

function. 
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PROBLEM SET 



1. a. The logistic population growth function satisfies the differential 
equation 

^ = Ax(N - x) 

where x-x(t) is the population size and N is the carrying capacity of the 
environment • Use this equation to show that tha maximum growth rate 
(max dx/dt) occurs at the inflection point x=N/2, Assume A is positive. 
(HINT: write the equation with v = dx/dt, and set dv/dx = 0). 
b. The ''relative growth rate" is defined by 

X ^ dt ^ 

where dx/dt is as given in part (a). Let x be given in "numbers of animals" 
and evaluate the units of R. The logistic curve is often used to describe 
"crowding effects" including intraspecies competition. Use the equation of 
part (a) to determine the value of x which maximizes R, and explain this 
result in terms of crowding, 

2, Leaves usually have small openings called stomata to allow passage of 
gases between their interior and exterior. Through them carbon dioxide 
passes in for capture by photosynthesis and the resulting oxygen passes 
out. Water vapor also escapes through the stomata, sometimes leading to 
dehydration. Thus plants have guard cells around the stomata to regulate 
their size. Action of the guard cells varies the shape of stomatal open- 
ings from a long narrow slit to nearly a circle. Throughout most of this 
variation the opening has approximately the shape of an ellipse with a 
constant length perimeter, typically about 35 y. 

24 
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2.(cont.) A good approximation to the perimeter of an ellipse is 

P = 27r/(a^+b^)/2 

The area is given by 

A = TTab 

Set P=35 and evaluate the area, A, in terms of just the width, b. Show 
that the area reaches a maximum when a=b, i.e., when the stomatal opening 
is a circle. 















k- 2a — 

1 


: -1 



Constant-perimeter ellipses. 



3. In the 0^ consiimption model represented by equation (7), the variable Q 

has units of yl/hr. Write Q in ml/hr and find for W=1000 mg. Wh^t are 

the units of -1^ .? What might this derivative represent biologically? That 
dm 

is, why would a biologist be interested in this derivative? 

4. A study of shape changes in nemerteans and flatworms (Alexander, 1968) 
theorizes that a basement membrane encloses the body and contains fibers 
which run in helices around the body. From Figure b, if the length D of 
the fiber is fixed, then £ = D'cosB and when the <5hape is cylindrical, the 
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4.(cont.) circumference and volume are 2Tir = Dsin6, v =•■ -nv^Z. 





2-nr 



a. Worm membrane 



b. "unrolled" 
membrane 
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a. Express v as a function of 6 (with D a parameter) and show v=0 when 
0=0 or e=7r/2 radians. What do these two cases mean physically? 

b. Find 6 which gives a maximum for v. Prove it is a relative maximum 
and not a relative minimum, N'ote that laboratory dissections show that in 
the relaxed worm the fibers run at about 55° to the axis of the body. Why 
would the relaxed worm have the maximum volume? 

5. Studies of insect flight (Alexander, 1968) use the theory of "forced 
vibrations" to explain the muscle action responding to nervous stimuli. 
If the "forcing function" is assumed to be 

Fsin(27Tnt), t = time, F = constant, 

then the steady amplitude A (magnitude of the vibrations) is given by 

A = F[(s-47r^n^m)2 + (27rnlO^]"'^ 

where K = viscous damping coefficient 
m = wing mass 
n = frequency 

s = stiffness of the vibrating medium 
Find the frequency (n) , called the resonant frequency, which gives the 

maximum amplitude. 

26 
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6. This exercise verifies the location of the relative minimum in the example 
of "water falling across a rock face". The function is 

z = x^ + - 3xy + 15 

Find all the critical points and determine which one is the relall-e minimu-n. 

7. Heat transfer in soils depends on many factors; among them is the variation 
in the soil itself (de Vries, 1975). 

Since most of the variation is in the vertical direction, a simple mathema- 
tical model is the one dimensional diffusion equation, 

^ 3t 3z 3z^ 

where we define 

T = temperature 
t = time 

z = vertical space coordinate 
C = volumetric heat capacity 
X = thermal conductivity 
When C and X are uniform in depth and constant in time, we have the 
simple diffusion equation: 

3T „ 3^T 
3t ^ 

where a = X/C is called the thermal diffusivity of the soil. The temperature 
at the surface gives boundary conditions for the model. For sinusoidal 
variation of surface temperature, we can write the boundary conditions as 

T(t,0) = T +6 costjJt 
a o 

T(t,a>) = T = constant 
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7.(cont.) Show that the solution to this model is given by the following 
function: 

T(t.z) = T + e ^^^^ cos(a)t - z/d) 
^ ' a o 



with 



d = (2a/Lo)^ 



Be sure to verify that this function satisfies the diffusion equation and 
the boundary conditions. 
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ANSWERS TO THE PROBLEM SET 



1. a. Find the x which gives max(dx/dt) . Differentiate dx/dt with respect 
to x» equate to zero and solve for x: 

^ # = A(N - 2x) 

0 = A(N - 2x) 
X = N/2 

Since the second derivative is negative, i.e. 

d^ (dx/dt) = -2A < 0 , 
dx^ / 

then at x=N/2, dx/dt is maximal. 

b. Find dR/dx, set equal to zero, solve for x: 

R = -^[AxCN-x)] = A(N-x) 

^ --^ -A < 0 
dx 

So no relative max exists, and the maximum must be at the lower end of the 
domain of x: since 0<x<N , then max (R) occurs at x=0. This model then 
implies that crowding effects are present whenever any animals exist. 

2. 35 = 27r/(a^+b^)/2 

.35.^ a"+ b " 

^ir 2 
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2. (cont.) The area A in terms of b is then 



Maximize A(b) : 



0 = - 4b^ , assume h ^ 0. 
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b = 4^ 



7T* 
2 



4b2 



2tt , since b > 0. 

Thus the maximum area occurs when a=b, i.e. a circle. 



3. To convert: (7) so Q is in ml/hr, we divide by 1000: 

Q = LSTxlO^^Cm+DW^ 



1^ = 1.87xl0"^W^ 
3m 

= .234 for W = 1000 

The units are then ml O2 per hour per unit of activity. 



4. a. First solve for r(0): i 

r = I;- sin e 



V = sine)^-(Dcose) 
v(e) « sin^ ecose 
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4. a. (cont.) 

v(0) =0 since sin(O) = 0 
v(90°) = 0 since cos(90°) = 0 

b. ^=7^ {2sinecos2e -sin^e} 

do HIT 

0 = 2sinecos^e -sin^S Assume e>0 
0 = 2cos^e -sin^e 
sin^e ^ 2cos^e 
tanS = 1/2 

e = 54.74° 

5, Maximize A by minimizing the denominator: 



0 = 4-[(s-4TT2n2in)2+(27TnK)2] = 2(s-47T Vm) (-87T2inn)+87T V 
dn 

0 = -2(ins-4TT2m2n2) + 
2ins-K2 = STT^m^n^ 
J _ 2ms-K^ 



6. First calculate the required partial derivatives:. 

f = 3x2 - 3y , I = 3y^ - 3x 



Now find the critical point (s) by simultaneously solving 

|S = 0 and 1^ = 0 
3x 3y 



Thus 

3x2-3y=0 3y^-3x=0 
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6.(cont.) From the first equation we obtain: 

Substitute into the second equation and solve for x: 

3(x^)^-3x = 0 
3x(x^-l) = 0 

Thus x=0, 1 . The minimum is said to be at (1,1) . 
With x=l we evaluate y: 

y = (1)' = 1 

To classify this critical point, we evaluate 



(6-l)(6-l)"(-3)^ = 36 " 9 > 0 

Since 



0 = 6 > 0 at (x,y) = (1,1) 



then (1,1) is indeed a relative minimum. 



7. First show the boundary conditions to be satisfied: 

T(t,0) = T -he e'^^^ cos(a}t-0/d) 
a o 

= T +0 coscot 
a o 

T(t,«>) = T + 0 e cos(a3t-«>/d) 
a o 

= T + 0 
a 

= T 

a 




28 



7. (cont'd.) Now show that the diffusion equation is satisfied: 



|r= -6 e'^'^'^o) sin(a)t-z/d) 
dt 0 



iJ. = _ 1 e ^-"^1^ cos(a)t-z/d) - e 3-)sin(a)t-z/d) 
dz d o o □ 



« -z/d 
0 e 



[sin(ti)t-z/d) - cos(ti)t-z/d)] 



1!t = (_ i) \^-z/d [sinCojt-z/d) - cos(iot-z/d)] 
dz d □ 



1 ^ 

(4) -^"^^ [cos(tot-z/d) + sin(tot-z/d)] 
d □ 



29 / , 

-r^'^"^ sin(iot-z/d) 
d 



Substituting into the diffusion equation: 

26' 

-e e'^'^Li sin(iot-z/d) == -a sin(iot-z/d) 

0 CI 

The equation is certainly true when sin(a)t-z/d) = 0. Now assume that 
sin(a)t-z/d)''^ 0 and divide both sides by -e^e"^^^ sin(a)t-z/d) : 

0) = 2a/d2 
Substituting for d, we obtain 



0) 



= 2a/C(2a/io)^] 



= to 



Thus the equation is indeed satisfied. 
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COMPUTER EXERCISES 

j 

Program DIFF is basically a plotting routine where the graph of a function 
demonstrates some property or use of derivatives. The program is designed to 
motivate the following concepts: the derivative as a slope, the zero slope 
at a critical point, partial derivatives, and the distinction between continuous 
and discrete rates of change. 

The specific user options and program features are detailed in the User's 
Guide for Program DIFF. The program is designed for easy conversion to be 
compatible with a continuous off-line plotter (such as Calcomp). Thus one 
option is setting XSLICE, YSLICE to "slice" and remove part of the three- 
dimensional graph to expose a hidden profile. This option is of limited 
benefit with a line printer, but is of great advantage with an off-line plotter 
for displaying critical points. 

In each exercise, the user must input certain parameter values- By 
repeating an exercise with different parameter values, the user can gain a 
be^tter intuitive understanding of how the behavior of functions and their 
derivatives depends on the chosen parameter values. 

1» In the example in the text treating oxygen consumption^ (see also problem 
3 above), the oxygen consumption (Q) depended on both the activity level (m) 
and the dry weight (W) : 

Q = .00187(m + 1)W*^ 

where the units of Q are mZO^/hr, For the computer exercise, we use the 
correspondence 

Q-^-z, m->x, W->y 
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so that the equation is 

z = . 00187 (x + l)y*^ 

Choose a value for y with 0 < y ^ 3000 mg and evaluate the first partial 
derivative 3z/3x at that value of y. Now use program DIFF: 

1. Choose function 1. 

2a. If a line printer is used, obtain the plot and check that the 

curvature with increasing x agrees with what your partial derivative 
suggests. Is the shape supposed to be a straight line, concave 
- upward, concave downward? 

2b. If an off-line printer is used, set XSLICE = 0 and YSLICE = your 

choice of y. Obtain the plot. Does the exposed profile agree with 
your partial derivative? Repeat the exercise with XSLICE = 0, 
YSLICE = 0. Does the graph have the same general shape as Fig. 7 
in the text? What is different? Now reverse the roles of x, y: 
pick X so that 0 £ x <^ 1, evaluate 3z/3y at that value, set 
YSLICE = 0 and XSLICE = your value for x and obtain your plot. 
Does the exposed profile agree with your function for 3z/3y? 
I That is, is the shape supposed to be a straight line, circle, 

parabola ... ? 

2. This exert ise illustrates how different parameter values can affect the 
properties of a function and its critical points. Function 2 is used in 
this exercise and is a simple polynomial in two independent variables: 

z = (Pl)x^ + (P2)y^ 

First calculate all the partial derivatives needed to locate and classify 
a critical point for general values of PI, P2. Now use program DIFF: 



1. Choose function 2. 
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2a. If a line printer is used, set values for PI, P2 (from -1.0 to 1.0) 
and obtain the plot. Locate and classify the critical point for 
these parameter values. Now check that the plot does show the same 
type of critical point at the same location given by your calculations. 
Repeat for several values of PI, P2 and note the dramatic change 
in type of critical point. How does the location change? 

If an off-line plotter is used, choose values for PI, P2, XSLICE 
and YSLICE and obtain the plot. Repeat the exercise with different 
values for XSLICE, YSLICE only, and obtain the plot. Classify the 
critical point as described in the text and compare w^^-h the plot. 
By repeating the exercise with different slicing values (XSLICE, 
YSLICE), you can search for the critical point by observing the 
changes in the exposed profile. For example, run the program four 
times using the following values: 

(XSLICE, YSLICE) = (-2 ,-2) (-2 ,-1) , (-2, 0) , (-2 ,1) 

You car. also keep YSLICE = -2 and vary XSLICE. Check that the 
exposed profile agrees with your first derivative. For example, 
if XSLICE = -2, YSLICE = -1, then the sliced edge should represent 
3z/9x at y = -1. 

3, 'When a beam of rectangular cross-section is cut lengthwise from a log, 
its strength can sometimes be well described by the following function: 



S = k W 



where 



S = strenpfh of beam 
W = width of beam 
D = depth of beam 

as shown in Figure A, and where, k is a constant which depends on the type 
of tree used. Assume that the log is circular in cross-section and that the 
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beam Is cut so each corner reaches the outside of the log, as shown. 
Assume k = 0.1. 




Figure A, Rectangular beam cut from a log. 

a. Find the width and depth of the beam which give the maximum strength. 
Assume the radius of the log is r. 

b. Now use program DIFF: 

1. Choose function 3. 

1 ^ ■ 

2. Let P2 = radius of the log. Choose P2 so that 0 < P2 < 10. 

3. Plot the function. . From the graph with y representing S and x 
representing W, estimate the width W which gives the maximum 
strength. 

4. Now calculate the depth D. Check that your estimates for W, 
D agree with your gpneral formula of part a. 

4. Animal populations newly introduced into a region have been observed to 
increase rapidly in number and soon thereafter to fall drastically (Caughley 
1970a, b) as indicated in Figure B. One theory is that the population at 
first "senses" an infinite food supply and then reproduces rapidly to the 
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point of overgrazing the area. The birth rate remains the same but the 
death rate (perhaps of young) dramatically increases until the population 
is low enough to match the food supply. The population then increases 
more slowly and seems to stabilize. Migration is also involved but is poorly 
understood. The second population rise to a "steady-state" suggests some 
adaptation or social "learning" by the population concerning their new 
environment. 



N 

N(o) y 

0 Time t 

Figure B. Population dynamics of introduced species. 

a. Consider the logistic equation written as follows: 

N(t) = K(l + be"^*^)"^ 

where K, b,' r are positive constants, and N(0) < K. In an attempt to 
describe the behavior in Figure B, one could modify the above logistic 
function by allowing r to vary over time, i.e. r = r(t). Can the population 
ever exceed K? If not, then this modified logistic function is not a good 
model. Answer this for a general r(t) in any fashion: by examining the 
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function N(t) itself,^ by evaluating the maximum of N(t), etc. 

An alternate' model is developed if the discreteness of the popula- 
tion is taken into account. The logistic function in part (a) assumes that the 
population size, N(t), will change continuously as time changes. In certain 
species, populations do not change size smoothly (see Sladen and Bang 1969), 
The breeding time is a short period, once a year so that many off-spring are 
born at the same time. The population size then increases in large jumps. 
This discrete change is not modeled well by differential equations. A closely 
related field is the study of "difference equations" (see Goldberg 1958) where 
discrete changes are allowed. The derivative of the logistic function satis- 
fies the differential equation 



- r»(t) 



K-N(t) 



K 



From this formulation, we see that the logistic model assumes that the 
population is always aware of how far from the carrying capacity is the current 
population size and so continuously adjusts for this difference. (Note that 
the derivati've depends on this ^^ference, K-N(t)). Actually, plentiful 
food one season may cause too many births the following season, with a definite 
time lag between abundant food and severe overpopulation and with few adjust- 
ments in between. 

The computer is used here to illustrate this inadequacy of differential 
calculus. The derivative above is replaced by a difference: 



N(t+1) - N(t) 

(t+1) - t 



fK-N(t)l 



K 

This model assumes the population changes only once per time period (e.g. , 
annually) and is unable to adjust in between. Is this more or less realistic 
than the original logistic model? 
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The computer program for this problem plots the population size (N) as a 
function of time (t) for a given set of parameters (constants). You should 
note in each plot whether N fluctuates at all, whether the initial drop is 
steep, and whether the population size eventually stabilizes at the carrying 
capacity (K) . 

r 

The input parameters are identified as: 
Pi = initial size 
P2 = growth rate (r) 
However, the computer program restricts PI to [-10,10]. The range of N 

o 

(initial population size) is chosen to be [50,250]. Thus we use PI to repre- 

sent, but not equal, the initial size as follows: 

N - 150 ' 

PI = ^ . 

10 

The carrying capacity is fixed at K = 500. Investigate the different 
kinds of behavior of the solution by using program DIFF: 

1. Choose function 4. 

2. Choose PI, P2 so that -10 ^ PI i 10 and 1 ^ P2 ^ 4. , 

3. Obtain the plot. 

Check that the initial size on the graph is correctly represented by your 
choice of PI. Repeat the exercise if you wish, but be sure two of your plots 
have PI = -10, P2 = 3.5 and PI = -10, P2 = 1.5. Can you estimate the value 
for P2 above which the oscillations are no longer regular (or periodic), 
i.e. May's "chaos'V(1978) ? 

Do you think this "difference" model could be a fairly good description 
of the introduced population? How would you include social or genetic 
"learning" into the model? That is, once the population size rises and then 
falls, what parameters might change in the model so that the next rise would 
gradual? 
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5. The reaction (Y) of the body to a dose (X) of drug can be represented by 
**he function: 



where PI and P2 depend on certain body characteristics and on the maximum dosage 
which can be administered. Y indicates the strength of the reaction, measured 
in millimeters of mercury if blood pressure is being tested, or perhaps degrees 
Celsius if change in body temperature is being measured. 

Find the dose that has "maximum sensitivity," i.e. where the rate of 
increase of Y is greatest. Is this "critical point" a maximum or minimum for 
Y*? yiiat is this point called on the graph for Y? 

Now use program DIFF: 

1. Choose function 5. 

2. Choose ri, P2 so that 0 ^ PI < 5, .1 < P2 < 1.0. 

3. Plot the function and calculate the correct value for X at the 
"critical point." Check that the graph of Y(X) has the behavior 
you predicted at this point. 

ANSWERS TO THE COMPUTER EXERCISES 

1. 32/3x » (.00187y'^) « constant for given y. 
Exposed profile should be a straight line. 



Y(^') « x2(Pl/2 - P2-X/3) 



PI 
2 




32/3y - (.00187) (x+l)(.7)y 



Exposed profile curves upward with decreasing slopcj. 
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2. 32/9X = 2(Pl)x, 82/3y - 2(P2)y 

32z/9x2 = 2(P1), 922/8x8y = 0, 32z/3y2 = 2(P2) 
32/3x = 0 if X = 0 
32/ay =0 if y = 0 . 

Thus (x,y) = (0,0) is the only critical point. 
To classify the critical point, evaluate: 

f32z 



32x 322 
3x^ 3y^ 



322 ]2 
3x3yJ 



4(P1)(P2) 



Pick PI = .6, P2 = 1.0 



Then A(P1)(P2) > 0, 



32z 



3x' 



= 2 (PI) > 0 and the point is a relative minimum. 




3. From the diagram, we use the Pythagorean theorem to obtain 

(2r)2 = w2 + d2 
Thus d2 = 4r2 - w2 

S = (.l)(4)r2w - (.1)W3 

= (.4)r2 - (.3)W2 

0 = (.4)r2 - (.3)W" 
The critical point is then at W = / 4r2/3 ' = 2r// 3 ', 
We have 



dis 

dW 



^ = - . 6W < 0 



so that when W = 2r// 3 , S is indeed at a maximum. 
The depth is then 



D = /4r2 - w2 
= / 4r2 - 4r2/3' 
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4. a. Even assuming r = r(t), N cannot exceed K if M(0) < K. We see this by 
writing N as a fraction. 
K 



N(t) 



1 +be-^^'^^'^ 



Since b > 0, the denominator exceeds 1, and thus 
N(t) < K 

b. Treat r as r(t) in the difference equation model, with 
r(0) = r^ 
r(l) » r /2 

0 

r(2) = r /3 

0 

r(3) = r(2) 

or some such scheme to decrease r as t±mi increases* 
5. Evaluate the derivative. 

Y' = (P1)X - (P2)X^ 
The "rate of increase" of Y is greatest when Y' is maximal. Find the 
inax(Y') by differentiating Y' and setting Y" « 0. 

Y«t ^ d(Y') ^ ^p^j _ (2)(P2)X 
dx 

0 = (PI) - (2)(P2)X 
X = (P1)/(2(P2)) 
The value for Y is then 

„ fPl l ^fPl (P1)(P2)1 ^ (PI) 3 
^ [2(P2)J [2 < ' 6(P2) J 12(P2)^ 

This point is an inflecti:T. point on the graph Y versus X. 
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USER'S GUIDE FOR PROGRAM DIFF 



Id entification 

DIFF - A program which displays properties of derivatives of mathematic 
functions 

Authors - Richard Hertzberg, Mark Bailey, Center for Quantitative Science 
in Forestry, Fisheries and Wildlife, University of Washington, 
Seattle. February 1979. 

P'rxrpose 

Program DIFF is the computer supplement to the instructional 
module "Calculus-Differentiation," by Richard Hertzberg, which reviews 
the basic principles and uses of differential calculus, with special 
emphasis on ecological and physiological applications. The computer 
prj:.;?ram displays graphs of seler^.ed functions so that certain derivatives 
are visible. Most plots serve as checks to the user's own calculations. 

Operation 

The user controls program DIFF through certain input variables. 
The principal one, NFC, selects the function to be displayed and also 
enables or disables other JLnputs listed in the INPUT TABLE below. Other 
inputs control the parameters in a function, the portions of a function 
to be displayed, the number of plots, and the structure of the plots 
themselves. 

Setting NFC « 1 selects the function 

F(x,y) = 0.n0187(x -f Dy^'"^ 
and disables the input variables Pi, P2, and NPLOTS. The user may alter 
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the variables XSLICE and YSLICE which act as "slicing" variables that cut 

through the 2 axis of the function aloag certain planes In order to reveal 

hidden profiles. The slicing effect merely sets the 2 values of the 

function to ZMIN (the smallest 2 value In the plot, which Is usually 

displayed as a blank) whenever x < XSLICE and y < YSLICE. 

Setting NFC = 2 selects the function 

J i 

2 2 
F(x,y) « p^x + p2y 

This option enables all of the Inputs listed In the INPUT TABLE. The 
parameters of the function, p^ and P2> are represented by the arrays PI 
and P2 which can hold up to six different sets of values and generate up 
to six different plots, where the number of plots Is controlled by NPLOTS. 
For example, If 

PI - 1, 2, 3, 4, 5, 6, 

P2 = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 

NPLOTS = 4 , ) 

I 

then DIFF would produce four plots as follows: 

plot 1: Ix^ + O.ly^ 

plot 2: 2x^ + 0.2y^ 

plot 3: 3x^ + 0.3y^ 

plot 4: 4x^ + 0.4y^ . 

Also, by repeating NFC = 2 with different values for XSLICE and YSLICE, 
the user can "search" for three-dimensional critical points by observing 
various function profiles. 
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Setting NFC » 3 selects the function 

- 0.1x(4p2 . Ji^ 
which can be written as 

F(x) « 0.l3p^(x) 

In terms of the notation used In computer exercise no. 3. This option 
disables the Input variables XSLICE, YSLICE, NPLOTS, PI and the last five 
elements of P2. P2(l) represents the function parameter 
Setting NFC « 4 selects the function 

^(y) = P2y(i-y/5oo) + y. 

This fimctlon represents the step-wise solution to the difference equation 

\+l " ^2^\S'^^^\i^^^^^ » where « y(Xj^) , 

which is the discrete form of the logistic differential equation 
dy/dx = P2y(l-y/500) , 

as given in computer exercise no. 4. This initial value for y CyCO)) 
is determined by the function 

y(0) = 10.0 p^ + 150.0 . 

This option disables the input variables XSLICE, YSLICE, NPLOTS, and 
the last five elements of PI and P2. Pl(l) represents the parameter p^^ 
in the initial value equation, and P2(l) represents the function parameter 
P2- 
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Setting NFC = 5 selects the function 



F(x) « (p /2)x^ - (p /3)x^ , 
1 2 



This option disables the input variables XSLICE, YSLICE, NPLOTS, and the 
last five c^lements of the arrays PI and P2, The function parameters 
and p^ are represented by Pl(l) and P2(l), 
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Program Organization 

The program fis organized according to the following flow chart: 



Start 



-> < ^ NODFLT = ? ^ ylgg- 



No 



Read default values for 
all parameters from the 
built in default file. 



Read in the user-supplied 
parameter values. 



FINIS = .T, 



Output the 
appropriate 
error messages. 



Yes 



No 



Check for errors in 
the user-supplied input 



Itiput errors found 



> 



No 



Calculate all the x,y,z points 
for the chosen functior*. All z 
values are set to ZMIN for all 
X < XSLICE, y < YSLICE 



Write out the x,y,z coordinates, 
for all points onto the plot file, 



Call the plotting subroutine 
PRNT3D which plots the output, 



Terminate 
program. 



Note : For 
NFC=2, this 
sequence is 
repeated NPLOTS. 
number of timer , 



EKLC 



is 
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Input 

All input is handled by the format free input package (Gales and 
Anderson, 1978) which permits a user to assign values to variables by a 
"name=value** convention. Not all variables need be explicitly assigned 
by the user, however, as unassigned variables automatically assume default 
values. The input consists of any number of data sets, each of which i^s 
terminated by a dollar sign ($) . Each data set generates a separate 
printer plot. 

The input for DIFF is divided into three classes: (a) variables 
having mathematical significance: NFC, XSLICE, YSLICE, PI, P2 and NPLOTS; 
(b) variables which control certain program operations, such as program 
termination or the handling of default input: IPRINT, ECHO, NODFLT, and 
FINIS; and (c) variables which control the printer plots (default values 
are in parentheses): X>IIN (0), XMAX (0), YMIN (0), YMAX (0), 2MIN (0), 
ZMAX (0), XRICH (0), YRICH (0), DFAULT (0), OVPRNT (.F.), AVE (.T.),' 
INT2D (.FJ, NX (60), NY (45), and 2MAP (0,1,2,3,4,5,6,7,8,9). The variables 
in the first two classes are explained in the following INPUT TABLE, whereas 
the printer plot variables are explained in the user's guide for PRNT3D 
(Gales > 1978). The user normally may ignore the PRNT3D variables since 
DIFF controls them internally. However, if he chooses to change any of 
them, he should do so with great care. In particular, the variables 
XRICH 3nd YRICH, if made too small, will cause DIFF to generate a verv 
large number of enrichment points, and consume far too much computer time. 



I 
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INPUT TABLE 



Name 


Type and 
Dimensions 


Range Limits 


Description 


NFC 


Integer 


1,5 


Identifies function to be 
plotted. 

Default value: NFC « 1 


PI 
P2 


Real (6) 
Real (6) 


-10,10 
-10,10 


Function parameters ^ Their 
physical significance and, 
to some extent, their range 
limits, depend on the 
particular function 
specified. 
Default values : 
PI = 1,1,1,1,1,1 
P2 = 1,1,1,1,1,1 


XSLICE 
Y SLICE 


Real 
Real 


XMIN.XMAX 
YMIN,YMAa 


Deletes part of the function 
Dy setting z.— z«riiiN lor 
X < XSLICE, y < YSLICE. 
Default values: 
XSLICE = 0 
YSLICE = 0 


NPLOTS 


Integer 


1,6 


The number of plots to be 
drawn. NPLOTS is used only 
when NrC^Z. 

Default value: NPLOTS = 1 


IPRINV 


Logical 


.F., .T. 


A logical value which 
causes the current values 
for all input variables 
(default as well as current 
user input) to be printed. 
Default value: 
IPRINT = .F. 


ECHO 


Logical 


,F., .T. 


A logical value which 
causes the user's input to 
be echoed if ECHO=.T., or 
suppresses echoing if 
ECHO=.F. 

Default value is: 
ECHO=.T. 
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Input Table (continued) 



Name 


Type and 
Dimensions 


Range Limits 


Description 


NODFLT 


Logical 


• F.J • T « 


A logical value which 
suppresses the input of 
default values if 
NODFLT=.T. 
Default value: 
NODFLT «= .F. 


FINIS 


Logical 


, F . , • T « 


A logical value which causes 
program termination if and 
only if FINIS=.T. 
Default value: 
FINIS=.F. 



Note : XMIN, XMAX, YMIN, YMAX determine the range of points to be 
plotted and are set internally for each function. 
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The last four variables deserve special mention • 

1. The logical variable IPRINT controls the output of all input variables 
which are currently in effect (default values as well as those specified 
in the current input set). Setting IPRINT=TRUE (or T or .T.) displays 
the input variables; setting IPRINT=FALSE (or F or .F.) sup^^resses the 
display. 

2. The logical variable ECHO controls the echoing of the input cards. 
Setting ECHO=TRUE causes the subsequent input set to the echoed; 
setting ECHO=FALSE suppresses the echo for the subsequent input set. 

3. The logical variable NODFLT can be used to inhibit the automatic 
assignment of default values to input variables. If NODFLT is set 
TRUE in the c urrent input set, then the current input set is assigned 
default values as usual, but all subsequent input sets merely 
accumulate more input values. In effect, the input values which exist 
after th i-th input set is read, become the default values for the 
(i+l)-th input set. The standard default values may then be restored 
by setting NODFLT=FALSE, but, again, the effects of this change are 
delayed until the next input set is read. To a limited extend, NODFLT 
permits a user to set up his own default values and can be very useful 
for executing a number of input sets which differ only in a few para- 
meters* Consider the following example in which a user wishes to 
slice the same function by using several different values for XSLICE 
and YSLICE; 

/INPUT SET 1: THE FOLLOWING VALUES BECOME THE DE FACTO/ 
/DEFAULTS FOR ALL SUBSEQUENT INPUT SETS: / 
NODFLT = TRUE, NFC = 2, PI = .1, .4, .5, 2.3, 

52 



-48- 



P2 = 2.2, 2.6, 2.5, 2.8. NPLOTS - 4, XSLICE = 0, YSLICE = 0, $ 
/INPUT SET 2: SLICE THE ABOVE FUNCTION/ 

XSLICE ■= -1, YSLICE = 0, $ 
/INPUT SET 3: SLICE IT ANOTHER WAY/ 

XSLICE ■= 0, YSLICE = -1, $ 
/INPUT SET 4: SLICE IT YET ANOTHER WAY/ 

XSLICE = 1.2, ySLICE = -1.2, $ 
/INPUT SET 5: NOW STOP/ 

FINIS « TRUE, $ 

A. The logical variable FINIS controls program termination. The user 
should add the card: 

FINIS = TRUE, $ 

as the very last input set. If , FINIS is not set, the program will 
terminate abnormally. 

Output 

DIFF produces sets of plots ^ via subroutine PRNT3D, which display 
function values. Each plot contains a title, legehd, x and y axis 
annotation, and printer plot lines or surfaces. The title displays some 
of the values, and ranges of values, for the variables used to generate 
the plot. 

The plot legend, in conjunction with the numbers along the x and y 
axes, allows users to interpret the plot numerically. The x and y axis 
numbers are of the form +N.NNN and differ from their true values by powers 
of 10 which are specified by the scale factors in the legend. For example, 
the, first line of the plot legend for RUN 3 reads. 
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SCALE FACTORS « X-AXIS :E+00 Y-AXIS :E-01 Z-AXIS:E+00 
hence the point, indicated by a "1" character near the top right of the plot 
is near .(x=l, 246, y=3,100) « (x=l. 246x10°, y«3,100xlO-^) « (x=1.246, y=0.31) 
The remaining two lines of the legend specify the number of points mapped 
to each z-level (-9 means > 99), 

Restrictions 

i 

The input values are restricted to the ranges given in the INPUT 
TABLE* Due to the physical interpretation of the function parameters, 
further restrictions are indicated in the description of each computer 
exercise. 

Error Messag es 

Three types of errors may occur in program execution: 

1. Syntax errors in the user's input. 

2. Parameter range exceeded. 

3. Plot parameter errors. 

Input errors 1 and 2 generate an appropriate error message, the calculations 
are skipped, and the next input set is read. For type 3 errors, the 
program suppresses the plot, outputs the error message and reads the next 
input set. If the plot file is empty, an error message is printed but the 
plot proceeds. Type 1 and type 3 error messages are listed in the user's 
guides for subroutines FFORM and PRNT3D. 

Sample Ru ns 

The control cards, input cards, and line printer output for five 
sample runs appear on the next few pages. 



0IFF#CH70000/T10/ P2 . 
ACC OUNT »♦•***♦ *» ♦♦*♦♦♦♦, 
COMMENT. 

C DM •^ENT. ♦*♦♦♦♦♦**♦♦♦♦*♦*♦•*♦*♦•♦*♦♦♦*♦*♦*♦<■****♦ ' 
COMMtNT.* THF ABOVE CARDS lOEKTIFY TH't JCR ♦ 
COMMENT.* (OIFF), SPrCIFY THE CFNTRAL ^F:^DkY ♦ 
COMMENT.* REOUIPEMENfS (70000 OCTAL)/. THF * 
COMMENT.* CENTRAL PROCESSUP TIME (10 5FC0N0S> * 
COMMENT.* THE JOB PRIORITY (P2)» AND THt * ' 

COMMENT.* ACCOUNT NUMBER AND PASSWORD ♦ 
COMMENT. *♦*♦♦♦♦♦♦♦♦♦♦**♦♦♦♦♦*♦♦♦♦♦*♦*♦♦♦♦*♦♦*♦♦♦ 
COMMENT. 

ATTACH^BDIFF/ID-QDIFF . 
ATTACH^ 8PR3D* ID-RPR3n . 

ATTACH, BFF,IO-?FF. 

COMMENT. 

COMMENT. ♦♦*♦♦♦♦♦♦♦♦♦'*'**♦♦♦♦♦♦♦♦♦*♦♦*♦♦*♦***♦♦♦♦* 

COMMENT.* THF Af^OVE CARDS ATTACH THE MAIN ♦ 
COMME.hT,* program CD IFF), ANT TWF SUPPORT' * 
COMMENT.* ROUTINES PRINT3D (BPR3D) AND TMf ♦ 
CO.rMENT.* FRFE FL'P", INPUT ROUTINE (3FF). TH£Y * 
COMMENT.* ARE ALL l\ ".INAf'Y ♦ 
COMMEHT. ♦*♦***♦**♦♦*♦*♦**'♦***♦*♦*• 
CDV;MENT. 

L0&D.90IFF, KPR3D, PF^. 
FXECUTc»DIFF . 
. CG.^.'-'FNT. 

CLEMENT. *♦*♦♦♦*♦♦♦♦'*♦♦♦♦■♦♦♦*♦♦♦♦♦♦**♦♦♦♦*♦♦*♦♦♦ 

CO^IMENT.* THE ABOVil CARDS LOAD THE ROUTINES ♦ 
COr'MENT.* INTO CrNT»AL MF^.O^Y AND PASS CONTROL * 
CCMMENT.* TO DIFF f O'i L J^fcCL'T I f;N * 
CG^I^E NT. ♦<■*♦♦♦*♦♦♦**♦♦♦♦♦**♦*♦♦♦♦♦♦♦ ♦♦♦♦♦♦♦♦*♦** 
CcrttENT. 
♦►r.R 

I 

I TML FOLLOWING Di^rAULT VALUES .APF USED 
/ 

/ NfC-l, NPLCTS-1» 

/ PI - 1.,1.>1.»1..1..I.» ^'2 - I. ,1. »!.>!. »].,!.. 

/ XSLICE - 0.> YSLICt ■ 0., 

/ ECrtG-.T., NODFLT-.F., FINIS*. F., 

/ Nx-bC, NY-A5f 2 SA ?-0, I, ?, 3, 4» , 6> 7, 0, 

/ xr.lN-O.O, XMAV-0.0* YMIN-0.0, YMAX-0.0. Z.'^IN-O.O* 7^Ay-0.0» 

/ ^AlCh-O.O, YRICH«0.0, DFAULT-0.0» 

/ OVf-SNT - ,S-,r AVE ■ .T.> INT2D ■ .f., IPPINT > .r. 

/ 

/ 

/ 

N F C 1 , S 

/ 

/ 

/44*t*«*«4i^*««ti*«*<=^*«««'«««««i4>PUN 2*************** '■♦***********■•*** 

/ 

NFC-2> S 

/ 
/ 

/«*«t44r«*««4'«*«44«t#«***««'«««««*RUN 3t««««*«**t-««<'«««*4«4*«««^#****** 
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' NFC-3, $ 

/ - I 

I A I 

/♦♦♦♦♦♦♦♦♦4*****.****************9UN C,***** ♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ **♦♦♦♦♦♦*♦*♦/ 

/ ' 

NFC"-^/ Pl--10# P?-1.2» I 
/ I 

I I 
****************** *****P.UH <5********************************/ 

/ f 
NFC-5, XPICH-O.0-j» Y^I CH-C . 00?.» IPRINT • .T., % 

I ' 
I I 
f*t*****»************************STCP**********************''*********^ 

I ' I 

PINIS-.T., J 

/ I 

I I 

/*4ft***********************************************************^***^*/ 

*f DP 
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OGRA.I -DIFF- READY FOR INPUT- 

THE: FOLLOWING DEFAULT VALUES AR£ USfD 



1.#1.>1>#1.>I>#1. 



NFC"1# NPLHTS-l* 
PI ■ 1.* l.» P? 
XSLICE - O.f YSLICE » 0., 
ECHO-.T.* NDOFLT-.F., FINIS*. r., 

NX«6C» 1* 2* 3»^* 5*ft» 7» 

Xf.IN«O.Oi XMAX-0.0# YMIN«0.0, YHAy«O.C, 2MIN«0.0> JMAy-O.O/ 
X3ICH»0.0» YRICH-0.0* DFAULT«0.0» 

OVPRNT « .F., AVE =" .T.» INT2D ■ .F., IPRINT • .F. 

NFC»1# S 
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XSLICE 
0 



OXYGEN CONSUMPTICIN 
FDR F « 0. 00187 * (I' + l.O) ♦ Y*»C.7 
0 YSLICE • 
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Y-AXISJ n+03 Z-AXIS: E*C? 

.''67(6t.)» .A00{-3), .f^3(77) 
.•V33( 3)^ 1 .067 ( 0 ), 1. 2"n( n ) 
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2** ♦♦♦♦♦*******♦♦♦♦*♦**♦♦ **♦♦****/ 



NFC-2/ S 
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PI Tr Pr AK, OP PASS CR ... 
FOR F ■ P1*X**.? ♦ P?«Y**? 
XSLICE • -2.0000 YSLICE • -2.0000 

PI - 1. 0000 P2 » 1 . 0000 
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V 

dec 



CALF FACTORS ■ X-AX1S« " + 00 Y-A^flSt t + 00 ?-AXISI E + O^^ 

0-ZA- .002(-9)> .^vl(7e)^ 1.779(79), 2.66e(77)f 3.5'^7(72) 

5-Z9- A.AA5(37)> ^. 33^(2?)* 6.2?3tli), 7.11K A)» P..noo( 1) 
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/ 



NFC-3# S 



ERIC 



PI 



STRFNGTH OP WOOD BEAM 
POH Y{X# D) - 0.1*X * D(X)**2 
l.OOOO P2 • 



1 .0000 



.229 



. 737 



.992 



1 .?A6 



1 . 50C 



3.100 



2.7^8 



2.39t) 



2.0A3 



1 .691 



1 .339 



.966 



.63A 



.282 



1 1 
i 1 
1 1 



OYl 

X 

0 



.2?v 



'l^l . 737 ' .992 I'Al'^ iJ.CC 



CALE FACTORS ■ X-AXI5« £+03 Y-AXISf E-01 
0-Z*.- 0(-9)» 1. 000 (20)^ ?,C00( C), 

t— 79- 5.C0a( 0), e.003 ( 0), 7.D00( 0), 



7-iXISi E+00 
3.C00( 0)# A.COOf 
b.OOOl 0), 9.000( 



0) 




62 



^'R 06 RAM -01 FF- READY FOP INPUT 

/4i4i4i4i4i4t4t«4t4t*4t*4E4t4c*4t*«4<°4i4'«>;'4i4i4><»4>*RUN 
NFC-'i* Pl«-10> P2-1.2# S 
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LC6ISTIC GROWTH - DIFFERENCE EQUATION 

FOP DY • ( P;?*Y* (1 .~Y/500 . ) )*DX 
PI - -10.0000 P2 - 1.2000 



7.000 



6.205 



5 .^39 



^ .61^ 



3.818 



P 

0 

P 3.023 

U 

L 

A 

T 

I 2.227 

N 



1 1.^32 

I 

E 



.636 



.305 
-X--" 



. 983 



1.3?? 



1 .661 

• — X — • 



2 .000 
X 



11 



1 1 



11 1111111111 11 111 in nil 111 1111111111111 



:3c: 



n^'i ."83 1.32? 

TIMF IN YFA5S 



•X--- 

,66 1 



SCALC FACTC9.S - X-A>ISi r+Oi Y-AXISl E+C-2 7-^^15 

20-Z^- 0(-9)# l.CO0(6F.), 2.0. 0{ 0), 3.000 0 , ^'S^^ 

2t)-Z9« ?.00C{ 0)/ 6.000 { 0)^ 7.000( 0)» B.000( 0)» 9.nD0( 



0 ) 
0 ) 
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^'?OGRAM -OIFF- READY FOR INPUT 

f ^^^^^49^^^ tttt ************ *****R[JH i*** ***************************** / 
NFC»5» XRICH-0.05. YRICH-O.OOt* IP^iNT • .T.# S 
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RENT VALUES FOR ALL INPUT VARIABLESi 



T5 



C= 
CF 

S 

LT 



5, 
If 

• lOOOOOOE't-Ol* 
.1000000E+01# 
.1000000F+01# 

• lOOOOOOE-^Ol* 



T# 



0 
0 



- • 0, 
4# 
8, 

0 , 

,1500000E+01# 
0 , 
■;2000000E + 00, 
0 , 
,9000000E+01, 
.,500C000E-01» 
.5000000E-02/ 




1^ 



op 



JSP- 



.1 OOOOOOE+01* 
.1000000E+01» 
.1 000000E*01# 
.lOOOOOOE+01* 



.1000000E-«-01# 
.1000000E-»-01# 



.1000000E-i-01# 
.1000000E^01# 



1# 

9, 



2# 
6> 



3> 
7# 
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BODY 
FOR Y(X) 
PI ■ 



RESPONSE TO FIXEOIORUG OOSF 
1 , GOOD ■ P2 • I .0000 



2.000Y 



1 .773 



1 .SIR 



1.091 



^6^ 



6 36 



^109 



162 



0 .229 .^63 .737 .992 
f X X X X— 



1 .?A6 
. — X- — 



1. 500 
--X 



11 111 



1 
1 
1 



1 
1 

1 



OYll 



1 1 



0 



.229 



.^,93 .737 .992 

DRUG DOSAGE 



1. 



X 

2A6 



•X 

5CC 



'^CAL'^ FACTORS - X-AXISt c+00 Y-AXISl E-01 
20-7^- 0{-9), 1.0C0(6O), 2.000( 0 , 

15-11" 5.0001 0), 6.003( 0)» 7.000( 0) , 



Z-AXIS I 
3.000( 0)f 

e.ooo( o)i 



E + GO 

4.00 0 ( 
9.000 ( 



0) 
0 ) 
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^/♦♦•♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦★♦♦♦♦♦♦♦♦STOP********************************'^ 
FINIS".!., $ 

P.A'^ -DIFF- TERMINATfO 
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